Observational Constraints on the Ultra-high Energy Cosmic Neutrino Flux 
from the Second Flight of the ANITA Experiment 
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The Antarctic Impulsive Transient Antenna (ANITA) completed its second Long Duration Balloon flight in 
lanuary 2009, with 31 days aloft (28.5 live days) over Antarctica. ANITA searches for impulsive coherent 
radio Cherenkov emission from 200 to 1200 MHz, arising from the Askaryan charge excess in ultra-high energy 
neutrino-induced cascades within Antarctic ice. This flight included significant improvements over the first 
flight in payload sensitivity, efficiency, and flight trajectory. Analysis of in-flight calibration pulses from surface 
and sub-surface locations verifies the expected sensitivity. In a blind analysis, we find 2 surviving events on a 
background, mostly anthropogenic, of 0.97 ±0.42 events. We set the strongest limit to date for lO'** — 10^' eV 
cosmic neutrinos, excluding several current cosmogenic neutrino models. 



I. INTRODUCTION 

The existence of cosmic -ray particles of energies above 
lO'^ eV, first established in the early 1960's, has become 
a problem of the first rank in particle astrophysics. Mod- 
els for their production must generate particle energies many 
orders of magnitude higher than achievable on Earth, and 
these models in turn require extreme source physics that has 
not yet been formulated in a self-consistent manner. Even 
more problematic, the propagation of such particles is limited 
by strong energy loss due to the Greisen-Zatsepin-Kuzmin 
(GZK) process fT]. Hadrons are produced via the Delta- 
photoproduction resonance by interactions with cosmic mi- 
crowave background (CMB) photons. This GZK cutoff in en- 
ergy limits the propagation distance of the ultra-high energy 
cosmic rays (UHECRs) to within 100-200 Mpc in the current 
epoch, and severely distorts the observed energy spectrum. 
Astronomy using charged-particle UHECRs is thus limited to 
the local universe, suffering from both the loss of source spec- 
tral information, and the difficulty in back-tracing UHECRs 
through intergalactic magnetic fields. 

At distances of several Gpc, corresponding to the star- 
formation maximum at redshift z ^ 1 , the higher energy and 
density of the CMB photons leads to even greater restric- 
tions on cosmic -ray propagation, and a more rapid energy 
loss to Delta photoproduction. However, information about 
the source particles does survive in the form of secondary 
neutrinos in the decay chain, known as the ultra-high energy 
(UHE) cosmogenic neutrinos, first described by Berezinsky 



& Zatsepin (BZ) Their momenta are unaffected by mag- 
netic fields, and they propagate without energy loss directly to 
Earth, retaining information about the cosmic distribution of 
UHECRs and their sources. 

The ANITA Long Duration Balloon experiment was de- 
signed to search for cosmogenic neutrinos via electromagnetic 
cascades initiated by the neutrinos in Antarctic ice, the most 
massive body of accessible, solid, radio-transparent dielectric 
material on Earth. We previously placed limits on the UHE 
cosmic neutrino flux from the first flight of ANITA (|3|], and 
provided a separate detailed description |4] of the ANITA in- 
strument, flight system, data acquisition, and analysis meth- 
ods. In this article we detail upgrades and augmentations be- 
yond the instrument and methodology previously reported. 



II. EXPERIMENTAL TECHNIQUE 

The second flight of ANITA (ANITA-E) launched from 
Williams Field, Antarctica on December 2 1 , 2008 and landed 
near Siple Dome after 31 days aloft, resulting in 28.5 live 
days. Fig. [T] shows an image of the payload on ascent after 
it had deployed to its full flight configuration, and an inset 
of the balloon and payload at float altitude. The mean ice 
depth in the field-of-view was 1 .4 km, approximately one at- 
tenuation length at sub-GHz radio frequencies |5]. ANITA-II 
flew at an altitude of 35-37 km above sea level (33-35 km 
above the ice surface), and was thus able to synoptically view 
a volume of ^ 1.6M km^ of ice. ANITA-II's sensitivity to 
cosmogenic neutrinos was improved substantially compared 
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to ANITA-I: the front-end system noise temperature was re- 
duced by 40 K, a 20% improvement in temperature [4]; 8 
additional quad-ridged horn antennas were added to the pre- 
vious total of 32; and the efficiency of the hardware trigger 
was optimized for impulsive signals. Also, the instrument 
was made much more robust to the effects of bursts of an- 
thropogenic radio-frequency (RF) interference with the abil- 
ity to mask channels from the trigger in the azimuthal sectors 
of the payload pointing at the noise source. Masking occurred 
on time scales of a minute. This upgrade significantly im- 
proved the livetime when in view of strong sources such as 
McMurdo and Amundsen-Scott Stations, and trigger thresh- 
olds remained at thermal-noise levels throughout the flight. 
The combined effect of all of these modifications led to an in- 
crease of about a factor of four improvement in the expected 
signal from typical cosmogenic neutrino models as compared 
to ANITA-I. 




FIG. 1 : The ANITA-II payload on ascent with the lower eight horn 
antennas deployed. The payload height is ~8 meters, and each an- 
tenna face is 0.95 m across. The inset shows the balloon and payload 
viewed telescopically at float altitude of 35 km. 

The ANITA-II hardware trigger selects impulsive radio sig- 
nals with broadband frequency content and temporal coher- 
ence over nanosecond timescales in the vertical polarization. 
The broadband nature of triggered signals was achieved by 
requiring power in multiple frequency bands, and temporal 
coherence was ensured by requiring impulsive power in a full- 
band channel. All simulations and laboratory measurements 
of Askaryan signals [6,-^ from ice sheet neutrino interactions 
show that the RF signal at the payload is predominantly ver- 
tically polarized, due to the surface Fresnel coefficient and 



Cherenkov geometry Htl- No strict pulse shape require- 
ments are enforced in the trigger to allow for variations in 
the shape of any individual neutrino-induced cascade, and the 
trigger is thus very inclusive. The trigger threshold rides at 
the ambient thermal noise level to maintain an approximately 
constant trigger rate of ^10 Hz, which is dominated (98.5%) 
by incoherent thermal noise fluctuations. These thermal-noise 
triggers have well-modeled statistical probabilities for pro- 
ducing random impulses, and are highly suppressed in analy- 
sis with requirements of spatial and temporal coherence. Most 
other triggers are from anthropogenic sources. Such signals 
can mimic neutrino-like impulses, and may arise from high- 
voltage discharges in electrical equipment or from metallic 
structures charged by blowing snow or related effects. To re- 
move such backgrounds, we identify active and prior human 
activity in Antarctica and optimize pointing resolution to reli- 
ably associate anthropogenic signals with known sources. We 
know of no expected particle-physics backgrounds. 

Extensive calibration and validation of the system response 
and trigger efficiency for a range of impulse signal-to-noise 
ratios (SNRs) was done both prior to flight and with ground- 
to-payload impulse generating antennas during flight. These 
in-flight pulser systems were located at the launch site and 
at a remote field station at Taylor Dome on the edge of the 
Antarctic plateau. In both cases, impulse generators were fed 
through an antenna immersed in the ice as deep as 90 m 101 • 
In-flight measurements of the impulses from Taylor Dome 
provided validation that refraction effects on signal propaga- 
tion through the ice surface do not significantly affect the co- 
herence of the received signal, to distances of 400 km. Mea- 
surements of trigger efficiency in-flight were consistent with 
expectations from the ground calibrations, considering the 
narrower frequency content of impulses from Taylor Dome. 
Ground-to-payload signals provided an equally critical func- 
tion in detector alignment and in determining the precision of 
directional reconstruction. ANITA's antenna signals are com- 
bined via pulse-phase interferometric methods fl, resulting 
in a radio map ("interferometric image") for each polariza- 
tion of the intensity as a function of payload elevation and az- 
imuth. The largest peak in either map determines the direction 
of the signal source. The payload coordinate frame is tied to 
the geodetic frame via onboard GPS and sun-sensor measure- 
ments to a precision of about 0.1°. The pointing resolution 
in our analysis is 0.2° — 0.4° in elevation and 0.5° — 1.1° in 
azimuth, depending on the SNR of the event. 



III. DATA ANALYSIS 

Our signal region comprises events that do not come from 
known sources of human activity ("camps") and have no 
geospatial partner event within distance and angular separa- 
tion criteria. We blinded the signal region in two ways si- 
multaneously. First, we inserted an undisclosed number (12) 
of neutrino-like calibration events at random, concealed times 
throughout the flight. Second, we blinded ourselves to events 
in the signal region i.e. events that have no spatial partner. The 
analysis efficiency and background are estimated before open- 
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ing this hidden signal box. Prior to unblinding, the complete 
set of in-flight calibration pulser data is used to optimize the 
pointing reconstruction of the analysis. In this section, we de- 
scribe the sequence of analysis cuts that we applied to the sam- 
ple until the final step, opening the blind signal box. Table U 
shows the the total event sample, the sequential number of sur- 
viving events as each analysis cut is applied, and the efficiency 
of each cut. The number of Hardware-Triggered Events is 
the total number of events recorded. Quality Events exclude 
on- and off-payload calibration pulser events, unbiased auto- 
trigger events, and a smaU fraction (^1%) of events with data 
corruption. 



TABLE I: Event totals vs. analysis cuts and estimated signal efficien- 
cies for ESS spectral shape 1 10]. 



Cut requirement 


Passed 


Efficiency 




Vpol Hpol 




Hardware-Triggered Events 


~ 26.7M 




(1) Quality Events 


~21.2M 


1.00 


(2) Reconstructed Events 


320,722 


0.96 


(3) Not Traverses and Aircraft 


314,358 


1.00 


(4) In Clusters <100 Events 


444 




(5) Isolated Singles 


7 4 


0.64 


(6) Not Misreconstructions 


5 3 


1.00 


(7) Not of Payload Origin 


2 3 


1.00 


Total Efficiency 




0.61 




FIG. 2: The vertical-polarization reconstruction position of all 
21.2 M Quality Events, with the color scale representing the aver- 
age peak value of the interferometric image for events which fall in 
that bin. The dark blue region is consistent with pure thermal noise. 
The green, yellow, and red regions are increasingly non-thermal and 
are consistent with anthropogenic noise. 



The great majority of remaining events are random thermal- 
noise-coincidence events in which individual waveforms are 
impulsive, but the ensemble of signals are incoherent and 



uncorrelated between antennas. We create an interferomet- 
ric image in each polarization by cross-correlating waveforms 
from neighboring antennas and summing the total normalized 
cross-correlation value for each elevation and azimuth. We 
construct a "coherently summed" waveform given the direc- 
tion of the largest peak in either map using the antennas that 
are closest to that peak. Fig.|2]shows the vertical-polarization 
reconstruction position of all 21.2 M Quality Events. The 
color of each bin represents the average value of the peak of 
the interferometric image for events which fall in that bin. A 
value of < 0.05 corresponds to the thermal-noise floor The 
clear separation of thermal and anthropogenic events moti- 
vates our accounting for each background separately. Nearly 
all of the non-thermal excess will be removed in analysis be- 
cause of its association with known bases. 

Thermal noise power originates from the ice and the elec- 
tronics. Because thermal noise does not correlate between an- 
tennas, the peak in the interferometric image is equaUy likely 
to be in any direction. We establish thermal noise cuts based 
on a sample of 2.1 M above-horizon (non-neutrino) events 
from a quiet period during the flight, when ANITA was not 
in view of McMurdo Station. To establish Reconstructed 
Events, we reject thermal noise by cutting on the peak value 
of the normalized cross-correlation in the interferometric im- 
age, the peak of the envelope of the coherently summed wave- 
form, and its fraction of linear polarization. We also require 
that events reconstruct to locations where there is ice on the 
ground, down to 35° below horizontal. Using a modest extrap- 
olation verified by simulation, we reject thermal noise events 
at a level of < 2.5 x 10^*^, corresponding to a background of 
0.50 ±0.23 events out of the 21.2 M Quality Events. The un- 
certainty represents a variation of fit parameters used for the 
extrapolation. For more detail, see Reference ifTTIl . Events 
which pass this set of cuts have been filtered for the pres- 
ence of continuous-wave (CW) interference - typically nar- 
row frequency-band signals from radio transmitters - and are 
largely impulsive in nature. 

Remaining events are mostly anthropogenic. To cut anthro- 
pogenic noise at a level of < 3 x 10^^, we require that any 
neutrino candidate be a single, geospatially-isolated event. It 
is extremely unlikely for two neutrino events to occur in close 
proximity to one another; correlated event locations are likely 
to indicate human activity. We first remove all events asso- 
ciated in both time and space with known flights and over- 
land traverses, so the remaining are Not Traverses and Air- 
craft. We cluster all remaining events with known active and 
inactive camps, other events that also pass all previous cuts, 
and locations where low-level RF power has been detected 
by statistical correlation of weak signals, shown in Fig. |2] 
Events are clustered with a base or another event if the dis- 
tance between the reconstructed location(s) is <40 km or if 
the angular separation between the locations is < 5.5 times 
the pointing resolution for the event(s). Table shows the 
number of clusters of reconstructed events vs. cluster multi- 
plicity for Clusters <100 Events to give a sense of the distri- 
bution of anthropogenic noise, and is used to estimate the an- 
thropogenic background contribution. Remaining unclustered 
events are Isolated Singles. The isolation requirement lowers 
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acceptance since each event removes a region of the ice sheet 
around it from the available neutrino target volume. We used 
seven largely independent methods to estimate that the an- 
thropogenic background remaining after our clustering cuts is 
0.65 ± 0.39 vertically -polarized (Vpol) events and 0.25 ± 0. 19 
horizontally polarized (Hpol) events ifTTll . (For example, in 
Table |II] there are 7 single events from known bases, 3 non- 
single clusters not from known bases, and 17 non-single clus- 
ters from known bases, yielding 3 x 7/17 = 1.2. After the 
polarization cut, this estimate becomes 0.91 in Vpol and 0.28 
in Hpol. Combining this estimate with six other estimates 
that use a similar technique gives the quoted value and un- 
certainty.) 

TABLE II: Event cluster multiplicities for all clusters with fewer than 
100 reconstructed events. All the cuts in Table |l] have been applied 
except for the Isolated Singles cut. The 5 events in the signal region 
correspond to the 2-1-3 events in the last entry in TableU 



Cluster Multiplicity Number of Clusters 



Camp Not Camp 



10-100 


8 


1 


5-9 


7 


1 


4 


1 





3 








2 


1 


1 


1 


7 


5 (Signal Region) 



At this stage, the polarization angle of the event is calcu- 
lated using the Stokes parameters, and the event is assigned 
to be Hpol (< 40°), Vpol (> 50°), or sideband (40° - 50°). 
There were no events remaining in the sideband. Events that 
are Not Misreconstructions have a low probability of mis- 
reconstruction; we remove any event that clearly peaked at a 
sidelobe of the pattern in the interferometric image. This cut 
was tested using events from the ground-to-payload calibra- 
tion pulsers and known camps. The requirement that events 
must be Not of Payload Origin removes events associated 
with local interference originating on the payload but missed 
in the quality event stage; these are easily identified. 

Upon opening the blind box, we first examine what hap- 
pened to the 12 undisclosed, inserted neutrino-like events. Of 
the 12 inserted events, 11 were unique events with one du- 
plicate event. Of the 1 1 unique events, 8 were Reconstructed 
Events, consistent with the calculated reconstruction analysis 
efficiency for such low-SNR events. Expected neutrino events 
have a wide range of SNR, but always at least as large as the 
inserted events. 



IV. RESULTS 

After all cuts are applied, two events remain in the Vpol 
channel, and three in the Hpol channel. After clustering cuts, 
the thermal noise background reduces to 0.32 ± 0.15 in each 
channel. The total background is 0.97 ± 0.42 events in the 



Vpol channel, and 0.67 ± 0.24 events in the Hpol channel. 
Thermal noise backgrounds are likely reducible in future anal- 
ysis. Fig.[3]shows the reconstructed locations of the remaining 
events (large blue squares) on the Antarctic continent and the 
payload position at the time of detection (small black square 
connected to the blue square), along with camps (red points), 
and locations of low-level RF noise (black points). 




FIG. 3: Events remaining after unblinding. The Vpol neutrino chan- 
nel contains two surviving events. Three candidate UHECR events 
remain in the Hpol channel. Ice depths are from BEDMAP L12il . 



All three Hpol events show characteristics which identify 
them as geo- synchrotron radio emission from UHECR air 
showers, reflecting from the ice surface (including sea ice), as 
described in our ANITA-I results L13]. While ANITA-I saw 
16 such events, the much smaller number of Hpol events seen 
in ANITA-II is due to the change of the trigger to favor Vpol 
events to maximize neutrino sensitivity. 

The two remaining Vpol events are of unknown origin. In 
Fig.|4]we show some of the characteristics of these events, in- 
cluding the waveforms, frequency spectra, and interferomet- 
ric images. Event 8381355 has been filtered using the adap- 
tive filter developed for the analysis and is highly impulsive, 
with a nearly flat radio spectrum. Event 16014510 shows a 
central impulse with some additional distributed signal within 
10-15 ns of the peak, and a frequency spectrum peaking near 
400 MHz. The event still passes aU cuts if the 400 MHz region 
is filtered by hand. The reconstructed directions are robust, 
supporting identification as isolated events. The waveforms 
and frequency spectra are within the range of simulated neu- 
trino events. Both events are consistent in their locations and 
amplitudes with distributions of Monte-Carlo-generated neu- 
trinos. We lack adequate statistics to identify these two events 
as a unique non-anthropogenic population. 

We proceed to set a limit including systematic errors ifTill . 
The largest systematic error is on the acceptance, which is cal- 
culated using two independent Monte Carlo simulations |4]. 
The two simulations typically differ by 20%, which we take as 
a systematic error. The uncertainty on analysis efficiency in- 
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FIG. 4: Top: Waveforms of incident field strengtfi for tfie two sur- 
viving Vpol events. Event 8381355 is shown filtered between 235- 
287 MHz to remove weak CW noise. Middle: Corresponding fre- 
quency power spectra with an average thermal-noise spectrum shown 
in blue for reference. Bottom: Corresponding interferometric images 
showing the pulse direction. The dashed line is the horizon. 



eludes two effeets: statistieal uneertainty on the efficiency cal- 
culation using calibration pulses, and a systematic error from 
the comparison of the efficiency on simulated neutrino events 
and calibration events. This uncertainty on analysis efficiency 
is 3%. The uncertainty on the background is discussed above. 
The inclusion of systematic errors only worsens the limit by 
about 10% because the procedure given in Ref. lfl4ll accounts 
for both signs of systematic fluctuations. 

Two other systematic sources are theoretical in nature. 
Changes in cross section, o, can affect the limit in two ways: 
increasing the neutrino-nucleon cross section increases the in- 
teraction rate, but lowers the solid angle due to Earth shield- 
ing (and vice-versa). The net result is an event rate which 
scales as a" '*^. For a choice of cross section different than 
used here ifTsll . our limit can be adjusted accordingly. Simi- 
larly, the event rate depends linearly on the Askaryan electric 
field f?*]. Including reasonable variations on these parameters 
would affect the limit by ^ 10%. 



V. DISCUSSION 

Our model-independent lfl6i [TtIi 90% CL limit on neutrino 
fluxes is based on the 28.5 day livetime, energy-dependent 
analysis efficiency (68%-42% from lO'^-lO^^^ eV), the aver- 
age acceptance from the two independent simulations, lill llsill 
and 0.97 ± 0.42 expected background events including the 
systematic effects described above. Relative to the revised 



ANITA-I hmit H [H] shown in Fig. |5l the expected Umit 
from this data, in the absence of signal, is a factor of four 
more sensitive. We set the actual limit, shown in Fig. |5] us- 
ing our 2 observed candidates. Because ANITA-II saw more 
than the expected background, the actual limit is only a factor 
of two better than ANITA-I even though the a priori sensi- 
tivity is four times higher for ANITA-II. The ANITA-II limit 
supercedes the ANITA-I limit and would not significantly be 
improved by combining the results. 
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FIG. 5: ANITA-II limit for 28.5 days livetime. The red curve is 
the expected limit before unblinding, based on seeing a number of 
candidates equal to the background estimate. The blue curve is the 
actual limit, based on the two surviving candidates. Other limits are 
from AMANDA |19], RICE |20], Auger (H, HiRes (H, and a 
revised limit from ANITA-I [18]. The BZ (GZK) neutrino model 
range is determined by a variety of models L10i , i23r - r28i1 . 



Table Unl gives integrated event totals for a range of cosmo- 
genic neutrino models with widely varying assumptions. We 
also include for reference the expected number of events for a 
pure power-law neutrino spectrum that matches the Waxman- 
Bahcall flux bounds for both evolved and standard UHECR 
sources Issll . ANITA-II's constraint on cosmogenic neu- 
trino models strongly excludes models with maximally ener- 
getic UHECR source spectra which saturate other available 
bounds |24, 26, 32]. These models generally assume very 
flat source energy spectra which may extend up to 10^^ eV; 
our results are incompatible with a combination of both of 
these features. ANITA-II is now probing several models with 
strong source evolution spectra that are plausible within cur- 
rent GZK source expectations ll24i l26l l28l [30l [3lll . some at 
> 90% confidence level. The ANITA-II 90% CL integral flux 
limit on a pure E^^ spectrum for 10^^ eV < isv < 10^^^ eV is 
E^Fy < 2 X 10"^ GeV cm"^ g-i 'j'jjese differential and 
integral limits, as well as the individual model limits above, 
are the strongest constraints to date on the cosmogenic UHE 
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neutrino flux. 
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TABLE IIL Expected numbers of events Nv from several cosmogenic 
neutrino models, and confidence level for exclusion by ANITA-II 
observations when appropriate. 



Model & references 



predicted CL,% 



Baseline models: 

Vai-ious riO, 23^ 24, 2a, 2^ 0.3-1.1 

Strong source evolution models: 

Aramo et al 2005 [26] 2.6 

Berezinsky 2005 [30] 5.4 

Kalashev et al. 2002 [24] 5.9 

Baiger, Huber, & Marfatia 2006 [2^] 3.6 

Yuksel & Kisder 2007 [31] 1.8 

Models that saturate all hounds: 

Ycshida a/. 1997 [32] 32 

Kalashev et al. 2002 [24] 20 

Aramo ef a/. 2005 [2&1 17 

Waxman-Bahcall fluxes: 

Waxman, Bahcall 1999, evolved sources [331 1.5 

Waxman, Bahcall 1999, standard [33J 0.52 



78 
97 
98 
89 



> 99.999 

> 99.999 
99.999 
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